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We consider an interesting and practically important case of elastic domain structure, which is 
the analogue of c/a domain pattern with 90° walls in perovskites, and is solvable analytically for 
arbitrary misfit strain. There is no critical thickness, below which the domain structure cannot 
exist, when the "extrinsic" misfit is zero and the domains are of equal width. At the boundary of 
polydomain-monodomain transition the period of the pattern diverges, as does the dynamic stiffness 
of the domain structure. It is unlikely, therefore, that one can achieve a softness of the dielectric 
response of the c/a elastic domains in ferroelectric-ferroelastic thin films. 
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Equilibrium domain structures in epitaxial ferroelec- 
tric films may appear even in the case of complete com- 
pensation of the depolarizing electric fields by the elec- 
trodes or finite conductivity of the film. This takes place 
if the film also behaves as a ferroelastic. Formation of a 
ferroelastic domain structure was first considered a while 
ago as a mechanism to relax the misfit imposed by a 
substrate ||l[. In widely studied perovskite ferroelectrics, 
which are improper ferroelastics, several elastic domain 
(and heterophase) structures were predicted One 
of the typical structures, so-called c/a/c/a... domain pat- 
tern, is of a special fundamental and practical interest. 
Firstly, the ratio of the widths of the c- and a-domains 
changes with external conditions (film thickness and mis- 
fit strain). At some parameters the domain pattern does 
not exist, i.e. there is a phase transition between multi- 
and monodomain states. Although the corresponding 
phase diagrams have been discussed by several authors 
(see [HQl and references therein), the behavior close to 
the phase transition was not studied. The model calcu- 
lations [|j have predicted the response of the c/a domain 
structure to be large, even giant for some specific val- 
ues of the ( "extrinsic" , see below) misfit w and the film 
thickness and some data was even interpreted in this 
way 1^. The enhancement of the response was expected 
near the border of existence of the c/a structure in the 
w — I plane, i.e. near the phase transition. Previous 
analyses, mainly numerical, have not actually specified 
the system parameters needed to realize the giant sus- 
ceptibility. We attempt, therefore, an exact solution of 
a simplest relevant model, analogous to the actual c/a 
structures. The present results put strict constraints on 
the system parameters where the softness of the dielec- 
tric response might be expected. It becomes clear that 
it would be very difficult to realize the domain structure 
with such a property. As a corollary, the existing data 



Q cannot be interpreted as being due to a large con- 
tribution of the elastic domain walls. We reveal in the 
present paper some unusual features of the polydomain- 
monodomain phase transition and show that the response 
of the system becomes more "rigid" in the vicinity of this 
transition. 

We present here the first analytical calculations of the 
energy of a domain pattern analogous to the standard 
c/a structure. The structure can be viewed, without loss 
of generality, as a result of a ferroelastic transition in an 
epitaxial film, which breaks the symmetry of the parent 
phase and the substrate. All the strains are considered 
in the reference frame of the high symmetry phase, and 
the z-axis is selected to be perpendicular to the plane of 
the film. We consider the system far from the transition, 
where the pattern consists of the domains having the 
spontaneous strain, m^^, of opposite signs, separated 
by the domain walls inclined at 45° with respect to the z- 
axis, as dictated by the elastic compatibility conditions 
and the lattice symmetry. Fig. 1. In ferroelectric per- 
ovskites this would correspond to the 90° walls separating 
the domains with the polarization parallel and perpendic- 
ular to the plane. The film consists of a— and c— type do- 
mains with the widths ai and a2, respectively. From the 
exact expressions for the free energy F [see Eqs.(||),(^], 
we calculate the period 2a — ai + 02 of the structure, the 
parameter of asymmetry i5 = (ai — 02) /2a, and the cor- 
responding stiffness k. The stiffness k is closely related 
to the "mechanical force constant" discussed in Ref. 
We shall consider two stiffnesses: the usually measured 
dynamic stiffness Koo, calculated at fixed period of the 
structure, and the static stiffness kq, corresponding to the 
situation where the period of domain structure is allowed 
to relax via creation or annihilation of the domain walls. 
The latter corresponds to very slow processes, making its 
observation very challenging. One certainly deals with 
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FIG. 1. Schematic of the ferroelastic phase transformation 
in epitaxial film, when the "extrinsic" misfit between the par- 
ent phase and the substrate is also present (top panel). As a 
result, the domain structure with the period 2a — ai + a2 is 
formed in the film of the thickness I (bottom panel). 

the dynamic stiffness Kqo in applications. Wc shall show 
that the two stiffnesses exhibit quite different behavior 
near the transition: the static one diminishes, while the 
dynamic one diverges. We also study the change of the 
stiffness with the film thickness and discuss the possibil- 
ities of obtaining the films with small dynamic stiffness. 

We consider first the case where the sole origin of the 
misfit is the spontaneous strain, u'^^ — u^^, appearing at 
the phase transition. We shall call this misfit "intrinsic" , 
to distinguish it from "extrinsic" misfit of any other ori- 
gin. The latter includes e.g. the misfit of parent phase 
with the substrate, and other, "noncritical" , strain com- 
ponents, which might appear during the phase transition. 
The effect of extrinsic strain is a focus of the present pa- 
per. Similar to previous authors |^ we neglect the un- 
compensated electric fields, which would only increase 
the stiffness. 

The method of calculating the energy of the domain 
structure is the same as in our previous paper ||^ but with 
different equations of state for the stress tensor compo- 
nents aik- 



= (A + 2^l) 



(1) 

(2) 



with axx and (Jyy components obtained from (^ by cyclic 
permutation of x,y,z. Here A,/i are the Lame coeffi- 



cients, and 



^ik 



is the elastic strain, with 



the components of the spontaneous strain. In contact 
with substrate, the film with symmetry-breaking misfit 
with the substrate must split into domains in such a way 
that the average in-plane strain will be zero, since the ho- 
mogeneous strain in the substrate would cost an infinite 
energy. Without the "extrinsic" strain the spontaneous 
strain would alternate from domain to domain (which 



all would have the same width) as u^^ = 



,0 — 



±Uq, 



yy 



= 0. On the other hand, with the "extrinsic" misfit 
strain w the distribution of the spontaneous strain in the 
domains would be 



w. 



(3) 



We assume for the substrate the same equations of state, 
as (|l|),(|), but with = = ul, = 0. 

To find the energy of the domain structure, we have 
to determine the energy of the homogeneous and inho- 
mogeneous (stray) stresses. The energy of the homo- 
geneous stresses {Fh) , which appears for domains with 
non-equal widths (i.e. for 5 0), can be easily found for 
an epitaxial film with any domain structure, provided 
that the elastic moduli in all domains are the same, and 
the non-linear effects can be neglected. We can readily 
find Fh from the general expression for the elastic energy 
Fei = -U dVa,kU°^ |,| as 



Fh 
A 
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= Ml 



XX yy 



= Mull [{S - Sof + w^{l - u^)/ul] 



(4) 
(5) 



where So ~ —w{l + i')/uq is the relative extrinsic mis- 
fit, V — A/2(A + /i) the Poisson's ratio, and M = 
2/1 (A 4- ^) / (A + 2/i) = £:/2(l - i/^), where E is the 
Young's modulus. The overbar here and below marks 
the averaging over the film (domain patterii). The con- 
dition Uzz — was used together with Eq.(|l|) to obtain 
for the planar stresses Uxx — '2M{exx + ^^vy)i ^yy — 
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^xx{x^ y), where Cxx = Uxx - u"^x = -''J'xx = -"^^oO - w, 
^yy ~ ~^yy ~ Eq.(^). We have also used the 

absence of the in-plane strains, Uxx = Uyy ~ 0, im- 
posed by the substrate, and an obvious relation u^^ = 
i(l -f S)(uq + w) + i(l - S){-uo + w) = uqS + w. Note 
that Eq. (Q) is rather general and can be applied to other 
kinds of domain patterns in the epitaxial films (cf . ) . 

The stray energy of inhomogeneous stresses in a stripe- 
like domain structure periodic in a;— direction. Fig. 1, is 
to be found from exact solutions for the strain field pro- 
duced by the domains. The pattern is defined by the dis- 
tribution of spontaneous strains u"^ (x) and (x) . The 
condition of the local equilibrium, d<7ik/dxk = 0, gives 
two sets of two equations for the displacement compo- 
nents Ux, Uz in the film and the substrate with the use 
of a standard relation 2uik = dui/dxk + duk/dxi. The 
boundary conditions are given by the absence of stresses, 
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Czz = cr^z = at the free surface of the film {z = I) and 
at z ^ —oo. Both the strain tensor and the displacement 
vector are continuous at the film-substrate interface. The 
original system of partial differential equations is reduced 
to a system of the ordinary differential equations with the 
use of the Fourier transform and then solved, as described 
in detail in Refs. [^[0], and we obtain the following simple 
expression for the stray energy per area A : 
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where k = nnl/a The series is calculated with the 
result 
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where Li„ (z) = E^i ^V*" = z$(z,n,l) b = 

2Til/a. For (5 = this formula gives the total elastic 
energy of a symmetric domain structure (a pattern of 
opposite domains of equal widths). We have found ear- 
lier a somewhat similar formula for another symmetric 
domain structure of ai\a2 type Let us consider the 
following general cases. 

Zero extrinsic misfit {w = Q) .— In this case the en- 
ergy of the homogeneous elastic field is simply Fh /A ~ 
MuqIS"^. The stray energy has a simple form for the fol- 
lowing two limiting cases. In the standard case of narrow 
domains (a <C /) and (5 <C 1 
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One has to add the surface energy of the domain walls 
to find the equilibrium domain width 



Fdn,/A = \/27Vc 



(9) 



where 7 = Mu§A/\/2 is the energy of the unit surface of 
the domain walls, and A is a characteristic microscopic 
length [0. The total free energy, F = Fh + Fstray^ Fdw, 
is minimal for symmetric domain structure (i5 = and 
OLi = 02 = a), with the standard (Kittel) domain width 

(cf. i) 
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Using this result, we see immediately that the high- 
frequency response of the domain pattern (at a fixed do- 
main width a) is 
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>0, 
(11) 



so the pattern is stable, although it softens because of 
the negative contribution to the stray energy from the 
terms oc (5^ in (^). Note also that Kqo increases with (5, 
and this, as we shall see, is a general result. 

If the elastic modulus /i were soft, the domains can 
become wide, a ^ I {b 1). There 
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and one finds a symmetric structure with the large do- 
main width H 



47re^/^/exp(7rA/8/), 



(13) 



which would be ^ Hf we had a substrate with A ^ dat 
(cf. the answer for a ferroelectric capacitor |]l^). This 
is contrary to a previously considered Kittel case with 
narrow domains ok <C ^. In the case of the wide domains 
the response softens considerably, but remains positive, 
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Koo = 2'kMuqI- > 0. 



(14) 



It is unlikely, however, that the softness of the present 
domain structure is of any practical importance. The 
problem is that even a small extrinsic misfit in very thin 
films will push the domain structure to the boundary of 
its existence , where it becomes rigid for experimen- 
tally accessible frequencies of external field, as discussed 
below. 

Non-zero extrinsic misfit {w ^ 0) .— The energy of ho- 
mogeneous stresses is given by Eq. (|]). To find the do- 
main structure close to the boundary of its stability we 
need all characteristics in the limit (5 — *■ 1, and a S> I 
{b <C 1), i.e. close to a monodomain state. The stray 
energy in this limit is 
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Note that Roytburd has approximated the numerically 
computed stray energy with the functional dependence 
(1 — (52) , i.e. without the important log term P,p^. 
Minimizing the total energy, Ftot ^ Fh + Fstray + Fdw, 
with respect to the half-period of the structure a, we 
obtain the equation 
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which has a solution 
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where we have omitted the terms cx In In {41 j tt A) ^ 
In (4^/7rA) . It shows that the domain period diverges as 
1/(1 — (5^) close to the phase boundary with the mon- 
odomain state. The total energy then takes the form 



Ft. 



Mu\lA 



{&-i^f + <i>{\~S^)-\{\-8'')\ (18) 



where from the equiUbrium value of b is to be found 
(subject to the constraint < 1), with the parameter 
(j) = (£A/27rO^^^ < 1, and £ = In (4//7rA) the logarithm 
of a large number ~ //A 3> 1. 

The transition into a monodomain state {6'^ — l) o*^" 
curs as a function of the extrinsic relative misfit at 
1^0 1 = 1 — 0- This condition suggests a critical thick- 
ness of the film Ic, where the polydomain- monodomain 
transition takes place at a given misfit So (i.e. the phase 
boundary in Sq — I plane), 
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(19) 



Thus, the domain structure may exist only at I > Ic (Sq) . 
This formula is obtained for Zc/A ^ 1, but it also cor- 
rectly gives an absence of the critical thickness of the film 
when the parameter Sq is zero {Ic = when Sq = 0). The 
absence of the critical film thickness has been suggested 
earlier for another domain pattern, ai|a2, from numerical 
computations 1^ . This is contrary to the speculations by 
Roytburd, who obtained Ic ^ Q for 6q = 0, apparently 
as an artifact of the employed approximations The 
critical point is approached linearly with 5q 



<5 = l-(l-0-Jo)M 



(20) 



We see that the slope is dS/dSo = l/cj) ^ 1, so that 
the approach to the critical point is very steep, it looks 
almost discontinuous. 

We readily obtain the high-frequency (measurable) dy- 
namical stiffness close to the phase boundary. 
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which diverges when 5"^ 1, i.e. close to the phase 
boundary. This is in striking disagreement with the re- 
sults of Pertsev et al. who claimed a softness of the 
c/a domain structure close to the phase boundary, which 
actually does not materialize. The static stiffness, which 
has no practical significance since it requires very long 
time to relax the pattern in order to optimize the number 
of domain walls, sharply vanishes with increasing relative 
misfit 



Mull 
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FIG. 2. The relative half-period of the structure a/aK, 
where aK, Eq.(|ic|), is the standard Kittel domain width, and 
the asymmetry parameter 5 = (oi — 02) /2a versus the rela- 
tive extrinsic misfit (top panel). Close to the boundary of 
existence of the domain structure, where a — > 00, the static 
stiffness kq vanishes, whereas the (measurable) dynamic stiff- 
ness «txj diverges (bottom panel). 



see the exact result in Fig. 2. Interestingly, close to the 
transition into monodomain state the system splits into 
two groups of wide (oi) and narrow (02) domains with 
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(23) 
(24) 



We see that the width of the wide domains diverges, the 
density of the domains walls decreases, but the width 
of the narrow domains, 02, remains small. The narrow 
domains are somewhat compressed compared to the stan- 
dard Kittel width or 

It becomes obvious from the present analysis that it is 
unlikely that one can succeed in making a "soft" domain 
pattern with the small m epitaxial thin ferroelectric 
films with elastic domains, suitable for the applications. 
In fact, the only way to do this would be to avoid any 
extrinsic misfit, i.e. to keep = 0, and to reduce the 
film thickness. However, in this case the interval of the 
extrinsic misfit strain, which allows the very existence of 
the domains, is very narrow. Then, the proximity to the 
phase boundary in systems with two kinds of inequiva- 
lent domains means that the stiffness rapidly increases, 
oppositely to what one actually desired. We have tacitly 
assumed above that the stiffness of the domain structure 
determines its dielectric response. Indeed, one can as- 
sume that there is a spontaneous polarization parallel to 
the film plane in one domain and perpendicular to it in 
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the other in the present c/a domain structure (as in c- 
domains in perovskites) . One can also assume, as in Ref. 

that S — Seq — IPsEext/ Koc in external field, where P, 
is the spontaneous polarization, Eg^t the external electric 
field, and Seq is the equilibrium value of S in zero field. 
We have shown earlier that this relation may not hold, 
and in some cases it strongly overestimates the dielec- 
tric response. Additionally, the neglected contribution 
of uncompensated electric fields would only make the re- 
sponse stiffcr. However, even this overestimate leads one 
to a conclusion that the dielectric response of the pat- 
tern with inequivalent elastic domains in the epitaxial 
thin films is actually suppressed. 
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